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We demonstrate cavity sideband cooling of a single collective motional mode of an atomic ensemble 
down to a mean phonon occupation number (n) m ia = 2.01q'3. Both (n) m in and the observed cooling 
rate are in good agreement with an optomechanical model. The cooling rate constant is proportional 
to the total photon scattering rate by the ensemble, demonstrating the cooperative character of the 
light-emission-induced cooling process. We deduce fundamental limits to cavity-cooling either the 
collective mode or, sympathetically, the single-atom degrees of freedom. 



Cavity cooling [THS] is unique among laser cooling tech- 
niques in that it is applicable, in principle, to arbitrary 
scatterers of light. The energy spectrum of the scattered 
field — which governs the cooling dynamics and equilib- 
rium temperature — is shaped by the cavity resonance 
rather than by the internal structure of the scatterer. 
Cavity cooling thus offers enticing prospective applica- 
tions, from preparing ultracold molecular gases [6] [7] to 
continuous cooling of qubit registers with far-detuned 
light [5]. In experiments to date [9"HPT]. cavity cooling 
of one atom [5] or ion [TU] is well described by a semi- 
classical model [2] [3] . In the case of an ensemble, the 
coupling of many particles to a single cavity mode can 
yield nontrivial collective dynamics [TTM14) . such as en- 
hanced cooling of the center-of-mass motion [11] . 

Ensemble cavity cooling (Fig. [I]) differs markedly from 
conventional laser cooling, where emission into a plethora 
of free-space field modes allows for simultaneous and in- 
dependent cooling of all atoms, or equivalently, all mo- 
tional degrees of freedom of the ensemble. In cavity cool- 
ing, a single collective motional mode C can be defined 
that is maximally coupled to the cavity [13J . while all 
other ensemble modes are decoupled from the cavity due 
to destructive interference in the light scattering from 
different atoms. The coupling of C to the cavity is co- 
operatively enhanced by constructive interference in pro- 
portion to atom number [H] [11] , allowing C to be cooled 
faster — and to lower temperatures — than a single atom. 

Pioneering experiments [131 114] have recently demon- 
strated that the cavity-coupled collective mode C can be 
studied using the concepts of optomechanics [15]. In- 
deed, the cooperative cooling of C — in the limit of weak 
mixing with other ensemble modes — is equivalent to the 
single-mode cooling [TBTfTT] of macroscopic mechanical 
oscillators [HI [19] by radiation pressure. Compared with 
solid-state mechanical oscillators, the collective atomic 
oscillator C inhabits a different parameter regime — of low 
mass and correspondingly large zero-point length — that 
may facilitate observing the quantization of mechanical 
energy [20]. Furthermore, the internal degrees of free- 
dom in an atomic ensemble constitute an extra tool for 
manipulating the motional quantum state. The collec- 
tive motion could, e.g., be squeezed by quantum state 




FIG. 1. (Color online) Ensemble cavity cooling. A probe 
laser (red) is placed at red detuning from cavity resonance to 
enhance anti-Stokes scattering into the cavity (blue), which 
cools a single collective mode C (solid green oval) at a coop- 
eratively enhanced rate 7 C . Single-particle modes can only be 
cooled by mixing (at rate 7 m ) with C. This differs from or- 
dinary laser cooling (inset), where free-space emission causes 
the atoms to be cooled independently. 

mapping from the ensemble spin [21] . 

To cavity-cool the single-particle degrees of freedom 
in the ensemble, mixing between C and other motional 
modes may be introduced by an anharmonic or inho- 
mogeneous trapping potential, or by collisions. While 
such cooling has been the subject of significant theoreti- 
cal studies, including detailed numerical modeling [H[S], 
experiments confirming the predictions are few. 

In this Letter, we cavity-cool and directly observe the 
relevant collective mode C of a trapped atomic ensem- 
ble. The rate constant of the cooling depends linearly on 
both photon scattering rate per atom and atom number, 
demonstrating that the cooling relies on the cooperative 
emission of light by the ensemble. Our results are well 
described by adapting an optomechanical model |16j to 
our system, where the mechanical oscillator C has a very 
small mass M = (10 -23 — 10~ 21 ) kg, a frequency of 500 
kHz (half the 1 MHz cavity linewidth), and a compara- 
tively low quality factor Q — 19. We verify the agree- 
ment with optomechanical theory for a wide range of 
collective- mode occupation numbers up to (n) ~ 10 3 , 
and we demonstrate cooling down to (n) m j n = 2.0^Q'g, 
close to the theoretical limit for our parameters. 

The optomechanical interaction Hamiltonian H in our 
system arises from a position-dependent dispersive cou- 
pling of the atoms to the cavity mode. Formally, H 
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describes the dipole coupling of N atoms with posi- 
tion operators Xi to light in a standing-wave cavity 
mode ("probe" mode, with wavenumber k and annihi- 
lation operator a) at large detuning A from atomic res- 
onance relative to the excited-state linewidth T. Adi- 
abatic elimination of the excited state yields H = 
hi}^2^ =1 s'm 2 (kxi)a^a, where ft = g 2 /A — with vacuum 
Rabi frequency 2g — represents the dispersive shift of the 
cavity resonance due to a single atom at an antinode, or 
equivalently, the ac Stark shift experienced by such an 
atom per intracavity photon. In our experiment, simi- 
lar to Ref. [T3], the atoms are trapped along the cavity 
axis in an optical lattice incommensurate with the probe 
mode. In the Lamb-Dicke regime, where the deviation 
xi = ii — £i of each atom from the local trap minimum 
at £i satisfies ((kxi) 2 ) <C 1, the Hamiltonian H can be 
written in terms of a single collective mode C of harmonic 
motion at the trap frequency ui t [13] , with position oper- 
ator X = N- 1 Y^Li sm(2k£i)xi [H]. In terms of X, 



H = hgXa'a, 



(1) 



where we have absorbed an overall shift 5un = 
fl sin 2 (fc£j) into the cavity resonance frequency. Eq. Ill 
represents the canonical optomechanical interaction |loF 
[T7] describing an intensity-dependent force of strength 
HQ = NMlk per photon, or equivalently, a cavity fre- 
quency shift QX proportional to X. 

For a probe laser detuned from the cavity line of width 
k, small shifts \QX\ < K yield proportional changes in 
intracavity and transmitted power. The A-dependent 
transmission can be used to monitor mode C, while the 
X-dependent changes in intracavity intensity — delayed 
by the cavity response — induce either cavity cooling or 
its reverse process, loosely termed cavity heating: specif- 
ically, the delay converts the position dependence into a 
velocity dependence of the force on the atoms, which ei- 
ther damps or coherently amplifies the collective motion 
depending on the sign of the laser-cavity detuning pi [3] . 

Viewed in the frequency domain, the dissipative pro- 
cess arises from unequal scattering rates on the Stokes 
and anti-Stokes sidebands due to the cavity resonance 
The full optomechanical Hamiltonian [T5], with the 
interaction term given by Eq. [I] predicts a cooling power 



P c = NT, 



cV E r (((n)\£ + \ 2 -((n) + l)\£_\ 2 ), (2) 



for a mean occupation number (n) of mode C; here 
r sc = (a i a)Tg 2 /A 2 is the photon scattering rate of a 
single atom at a probe antinode into free space, rj — 
4<? 2 /(«r) the cavity-to- free-space scattering ratio (single- 
atom cooperativity) [3], E T = h 2 k 2 /(2m) the recoil en- 
ergy for atomic mass m, £ = TV -1 ^ sin 2 (2fc^), and 
C^. 1 = 1 =F 2i(S ± 0Jt)/n, where 8 is the probe-cavity de- 
tuning. In our experiments, where the atomic cloud is 
long (w 1 mm) compared to the 5-//m beat length be- 
tween trap and probe light, £ = 1/2. For u> t > k/2, the 
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FIG. 2. (Color online) Mean occupation number (n) of mode 
C vs. time during cavity cooling at r sc = 1.1 x 10 5 s _1 (gray 



squares), 2.3 x 10 J s _1 ( green circles), 



(gold 



triangles), and 6.4 x 10 J s~ (red diamonds). Each dataset is 
obtained by averaging variances from 10 traces. Inset: single 
trace of cavity transmission during cavity heating (t < 0, blue 
background) followed by cooling (t > 0, red background). 



cooling rate is maximized by placing the anti-Stokes side- 
band on resonance, 5 = —u>f Eq. [^indicates a collective 
rate constant 7 C = dP c /d((n)huj t ) that is proportional to 
N due to cooperative scattering: the larger the ensemble, 
the faster C is cooled. 

We study the cooling in a symmetric near-confocal op- 
tical cavity with linewidth k = 2tt x 1.01(3) MHz at the 
wavelength 27r/fc = 780 nm of the 87 Rb D2 line, mode 
waist w = 56.9(4) ^m, and cooperativity rj = 0.203(7). 
We trap 10 2 -10 4 atoms of 87 Rb in the state \5 2 S 1/2 ,F = 
2, nip = 2) in the cavity mode in a standing wave of 851- 
nm light, with trap frequency w t /(27r) = 480(40) kHz 
and typical trap depth U /h = 18(3) MHz. A de- 
polarized 780-nm probe laser drives the cavity on a 
TEM 00 mode at a detuning A/(2tt) > 70 MHz from 
the |5 2 S' 1/ 2,F = 2) -> \5 2 P 3/2 ,F' = 3) transition with 
linewidth T — 2tt x 6.1 MHz. The atom number N 
is measured via the average cavity shift Sljn ^Z\. To 
perform cavity cooling/heating, we detune the laser by 
5 = T K /2 ~ T^t from cavity resonance, simultaneously 
probing the position X via the transmitted light. Note 
that we work with blue light-atom detuning A > 0, where 
free-space scattering results in Doppler heating. 

We first verify cavity heating of mode C by choosing 
the probe-cavity detuning 5 = +n/2. Suddenly turning 
on the probe light triggers a collective oscillation that is 
rapidly amplified by parametric instability (inset to Fig. 
[2]). After typically 10 fis of this heating, we switch to 
cavity cooling at i5 = —k/2. The mean occupation num- 
ber (n) of C is obtained from the observed time trace of 
the transmitted photon rate R via the fractional variance 
a 2 = R 2 / R — 1 in a sliding 2-/is window. The linear ap- 



proximation X oc R — R gives the relation 0—0 
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FIG. 3. Collective cooling rates, (a) Energy decay rate 7 exp 
vs. scattering rate T sc for: iV=8000(700), A/(2tt) = 270 
MHz (solid black squares); iV=2800(400), A/(2tt) = 140 MHz 
(open red circles); iV=700(200), A/(2tt) = 70 MHz (solid blue 
triangles). Lines are fits to data, (b) TV-dependence of cooling 
rate normalized to single-atom scattering rate into cavity. 
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FIG. 4. Spectra of fractional transmission fluctuations Si /I 
taken with (a) N=230(50) atoms and (b) N=450(90) atoms 
during cavity cooling (S = — k/2; red lines) or heating 
(S — +k/2; blue lines). Black curves are fits; subtraction 
of the background level Sbg (gray curves) yields collective 
mode occupation (n) at S = ±k/2: (a) {n) + = 4.4 ± 0.7, 
(n>- = 2.3lg;5; 0>) (n> + = 7± 1, (n)" = 2.0l°;§. 



8(GX /k) 2 \C+ - £-\ 2 (n), where X = y/h(/(2Nmu! t ) 
and a\ is a constant technical- noise offset !22j. Fig. 
[2] shows (n) vs. time at four different probe powers, 
with fixed atom number TV = 2800(400) and detuning 
A/(27r) = 140 MHz from atomic resonance. The cooling 
is well described by an exponential decay with rate con- 
stant 7cxp that depends on the probe power. Consistent 
values of 7 cxp are obtained by fitting an exponentially 
decaying sinusoid to the averaged transmission trace. 

To compare 7 OX p to the predicted cooling rate con- 
stant 7c, we measure the dependence of 7 cxp on the pho- 
ton scattering rate r sc = i??7r 2 /(2A 2 ) per atom into 
free space for various probe-atom detunings A and atom 
numbers N. As Fig. |3ja) shows, the data are consis- 
tent with a linear model 7 exp = f(N)rjT sc + 7, n . The 
offset 7 m = 1.6(6) x 10 5 /s indicates a quality factor 
Q = iOt/lm ~ 19 for mode C, largely attributable to 
mixing with other motional modes in the anharmonic 
trapping potential. Note that our system allows cavity 
cooling at very low Q compared to solid mechanical oscil- 
lators (THJ |T9| because the "thermal bath" comprising the 
other 3N—1 ensemble modes has a sub-mK temperature. 

To verify the cooperative nature of the cavity cool- 
ing of mode C, we plot in Fig. (3^b) the fitted slopes 
f(N) as displayed in|3ja) vs. atom number N. Account- 
ing for the slight (< 20%) reduction of the cooperativity 
r\ due to atomic absorption, the measured dependence 
dj c ^p/d(i]T sc ) — 3.4(5) x 10 _3 A agrees well with the pre- 
diction from cavity cooling 7 c /(?7r sc ) = 3.0(2) x 10~ 3 A. 
This confirms that the collective-mode cooling speed in- 
creases linearly with ensemble size and is proportional to 
the total power scattered by the ensemble into the cavity. 

To determine the equilibrium temperature of C under 
cooling, we require — given our detection noise — a longer 
observation time than shown in Fig. [2] We therefore ob- 



serve the cooling or (for comparison) heating in spectra 
obtained from 150 time traces of the cavity transmis- 
sion, each 440-/JS long, with R = 1.2(2) x 10 9 s -1 . Fig. 

[4] shows normalized one-sided spectral densities Si /I 2 
of photocurrcnt I oc R with (a) N = 230(50) and (b) 
N = 450(90) atoms at a detuning A/(2tt) = 70 MHz 
from atomic resonance. Each spectrum displays a peak 
at io t with an area approximately proportional to both 
atom number N and mean occupation number (n) . The 
disparity in area between cooling and heating increases 
with N due to the cooperative nature of the processes. 

We fit the spectra in Fig. [4] with a quantum mechan- 
ical model (black curves) adapted [22] from Ref. [IB] . 
2 

The model Si /I = S mec h + S^g contains the signal 
•Smcch ~ (2^/k) 2 \C+ — £-\ Sx arising from atomic mo- 
tion with spectral density Sx] an d a background Sb g 
(gray curves) that is dominated by electronic photode- 
tector noise but also accounts for photon shot noise, 
slightly smaller fluctuations from laser phase noise, and 
frequency-dependent correlations between light noise and 
atomic motion. These last are responsible for the dips in 
Sbg below the white noise [19] . With the photon rate R 
and optomechanical coupling Q constrained to their in- 
dependently measured and calculated values, the cooling 
spectra are well fit by taking the collective mode to be 
coupled to a white Markovian bath with (ribath) = 3.1(4) 
motional quanta per mode; the corresponding coupling 
rate 7^ = 2.6(1.1) x 10 5 s" 1 is consistent with the mix- 
ing rate 7 m from Fig. [3] Fits to the heating spectra, 
complicated by sympathetic heating of other modes, in- 
dicate a higher mixing rate of 4.8(5) x 10 5 /s. 

The bath occupation is consistent with a measured 
upper bound on the axial temperature of 150(50) /iK, 
corresponding to (ribath) = 6(2) [22] . The white spec- 
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trum of nbath is a simplistic ansatz but helps to estab- 
lish the background Sbg and thus the motional spectrum 
Sx- By subtracting Sb g from the measured spectrum, 
we obtain a minimum mean occupation number of C of 



= 2.0 



+0.9 
-0.3 



with N = 450(90) atoms. Note that 
failing to account for the dip in Sbg would underestimate 

(^)min • 

We now consider limits to cooling the collective mode. 
For 7 C 3> 7m, the cooling power P c oc Nrj competes only 
with the iV-independent recoil heating P rcc w E T T SC of C, 
yielding a fundamental limit (n) > n + D(l + n ) / (Nrj) , 
where no = {n/Auit) 2 and D is a prefactor of order unity 
[3] . Thus, for large collective cooperativity N-q ^> 1 (easy 
to achieve), the resolved-sideband regime no < 1 in prin- 
ciple allows ground-state cooling [31 [THl E] of C. The 
thermal heat load from other modes mixing at rate 7 m 
with C then sets the limit (n) > (nbath)7m/(7m + 7c)- 
While this limit improves with increasing cooling rate 
7c for the values (tjj, k,Q — i^t/Ym) m Fig- @ amplifica- 
tion of low-frequency noise on approaching the regime of 
static bistability 7 C > uj^/k jTB] sets a bound (n) > 1.5, 
even though n = 0.3. 

A low occupation (n) of C is disadvantageous for cool- 
ing the individual atoms, since the absolute cooling power 
is proportional to (n) (see Eq. [2]). Cooling of all degrees 
of freedom is thus facilitated by strong mixing 7 m 7 C 
that keeps C in thermal equilibrium with the other 37V— 1 
modes. The cooling power per atom P c /N then ap- 
proaches that of an isolated atom. Thus, even in an 
ensemble, recoil heating sets a limit for the tempera- 
ture of individual atoms (n.j) > ?io + 1/rj that depends 
on the single-atom cooperativity rj: ground-state cooling 
requires rj > 1 [3 . Whether the same result holds in 
other cooling geometries, e.g., with transverse pumping 
13 [23] , is under investigation [23] . 

Even for rj < 1, ground-state cooling of C alone — 
in future experiments deeper in the resolved sideband 
regime — may enable the preparation of non-classical mo- 
tional states [HI [21] • Further, the sensitive detection 
demonstrated here for X can alternatively be applied to 
measure X 2 and thereby observe phonon shot noise [24] 
or perhaps even quantum jumps in n [20] , 
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In this supplement, we begin in Sec. I by reviewing the description — previously derived by Murch et al. in Ref. 
[1] — of a trapped atomic ensemble as a single-mode mechanical oscillator coupled to the cavity. In Sees. IIA-IID, 
we calculate the spectrum of intensity fluctuations at the output of the symmetric optical cavity, and we relate this 
spectrum to the spectrum of position fluctuations of the mechanical oscillator itself. We closely follow the approach of 
Marquardt et al. [2] for calculating the displacement spectrum of the mechanical oscillator, and wc additionally include 
lowest-order effects of laser phase noise [3]. After accounting in Sec. HE for technical effects in our photodetection of 
the cavity transmission, we derive in Sec. II F the relation of measured photocurrent fluctuations to the occupation 
(n) of the collective mode. Our calibration of atom number is described in Sec. III. Finally, in Sec. IV we describe a 
measurement of the thermodynamic temperature associated with all N modes of the ensemble's axial motion. 

I. HAMILTONIAN AND COLLECTIVE MODE 

We first consider the Hamiltonian H sys of an ensemble of N atoms that are harmonically trapped, with identical 
trap frequencies u>t, at various positions along the cavity axis. For dispersive coupling of the atoms to a cavity 
mode ("probe" mode) of frequency cj with annihilation operator a, we have 

+ huiQa) a. (1) 

Here, the atom-probe interaction is quantified by il = g 2 /A, in terms of the vacuum Rabi frequency 2g and detuning 
A > T of the probe mode from the atomic resonance with lincwidth T. In the Lamb-Dicke regime, where the 
motion x; t = ii — £j <C k^ 1 of each atom in its trap is small compared to the probe wavelength, the optomechanical 
interaction term can be expressed in terms of a single position variable X = N^ 1 J2iLi sin(2fc^)xi corresponding to 
the collective mode C discussed in the main text. The momentum conjugate to X is expressed in terms of the single- 
atom momenta pi as P = N J2iLi sm (2fc£j)pj/ X)i=i sin 2 (2fc£i), such that X and P obey the canonical commutation 
relation [X, P] — ifi. To describe the full motion of the A-atom ensemble, one can construct an orthogonal basis 
comprising X, P, and additional coordinate pairs representing 3A — 1 other modes of ensemble motion with energy 
H±; in this basis, the system Hamiltonian becomes 

H sys = \^Mu? t X 2 + f- + h(uj + 8u; N + QX} a* a + H ± , (2) 

where M — mN 2 / J2iLi sin 2 (2fc^) represents the effective mass of the collective mode, Sujjy = f2 J2iLi sin 2 (fc^) is 
an overall shift of the cavity resonance due to the atoms, and Q — N£lk is the change in this cavity shift per unit 
displacement of the collective mode. 

In Eq. 2, the collective mode is entirely decoupled from all other ensemble modes due to our approximation of 
perfectly harmonic and homogeneous trapping. We allow for corrections to this simplified model by introducing a 
term H lm representing coupling of the collective mode to a bath, which, in addition to including effects of mixing 
with other axial modes, might include effects of radial motion or light-induced effects beyond the optomechanical 
interaction included in H sys . Further accounting for a coupling H K of the intracavity field to input field modes with 
energy i/drivc, the full optomechanical Hamiltonian takes the form 

H tot — h [uj + Suj n + 0X0(0* + c)] (a' a - (a^a)) + tko t c'c + H dlivc + H K + H lm , (3) 

where we have subtracted an offset associated with the average intracavity light level and absorbed the associ ated 
force on the atoms into a redefinition of the trap centers & . In terms of the zero-point length X a = y/Ti/ (2Mu ( ), the 
annihilation operator c for mode C has the usual definition such that X = X (c < + c) and P = iMuj t X n (c < — c). 
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II. MECHANICAL MOTION AND TRANSMISSION FLUCTUATIONS 

A. Equations of Motion 

Equation 3 is the standard optomechanical Hamiltonian [2, 4] describing a cavity mode, with operator a, whose 
frequency is shifted in proportion to the position Xq($ + c) of a mechanical oscillator of frequency w t . The last three 
terms in Eq. 3 can be written out explicitly using the standard input-output formalism of quantum optics [5, 6]. One 
thereby derives equations of motion for a and c in terms of the optical input operators a ln j, corresponding to the 
fields driving the cavity from its two ends labeled by je 1, 2, and a mechanical input operator c; n corresponding to 
the thermal bath. We shall describe the evolution of the field operators in a rotating frame at the drive frequency ujl, 
writing the intracavity field operator a = (a + d)e~ l0JLt in terms of a c-number a and a noise operator d that accounts 
for small deviations from the classical value, and similarly letting d m j = (a,i n j + d- m j)e~ %WLt . The classical field values 
are then in a steady-state relation 

= (iS - ^)a - ^ (a in i + a in2 ) , (4) 

while the deviations obey the linearized equations of motion [2] 

d = (i5 - ^)d + ia(c + c f ) - ^ (d inl + d in2 ) , (5) 

c= (-iw t - ^y)c - y/^c in + i (a*d + aS^j , (6) 

where S = ujl — (uj + 5ljn) represents the detuning of the drive field from cavity resonance for X = 0; and a = —QX n a. 

Defining the cavity response function XkC 1 ^) = — i(oj + 5)] and the mechanical response function Xm(^) = 

l/[-f m /2 — i(u — u>t)], and rewriting the equations of motion in the Fourier domain, we have: 

J( -<„)]- 



d/x K (u) = ia[c{u) + c f (-w)] - J ~[d inl (uj) + d in2 (oj)}, (7) 
c/Xm(w) = i[a*d(oj) + a<?(-w)] - y/j^c in (uj). (8) 

Here, O(co) = f^^/v ^^Oit) dt/y/T denotes the windowed Fourier transform [6] of an operator O, and in calculating 

spectra (o(u))0(— u)j (for Hcrmitian O) we shall always implicitly take the limit T — > oo of a long window. 

We can solve Eqs. 7-8 to obtain the dependence of X(lo)/Xq = c(u>) + c)(—u) on the optical input fluctuations 
d m = (dini + ^in2)/v / 2 and the mechanical bath operator c- ln : 



-VT^HXm^-^CinM +Xm 1 ( w )Ci n (- W )] ~ ^V^^t a* \k (uj)d in (w) + «Xk (~ v)dj n (- U}) 

c» + *(-«) = ^ L (9) 

where 

AT(u>) = Xm 1 (^)x m 1 *(-^) - 2t \a\ 2 utUiu) (10) 

and 

n(w) = - x:(-w)- (ii) 

B. Input Field 

We shall allow one end of the cavity to be driven by a laser at frequency ljl that may have some phase noise. To 
account for this, we let <ii n i = <ii n o + ij3a- m i, where <ii n o represents quantum fluctuations, while 0(t) <C 1 is a real- 
valued stochastic variable representing the phase noise. At the other end of the cavity, we will admit only vacuum 
fluctuations (setting a; n2 = 0). 
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The quantum fluctuations (j = 0, 2) satisfy 

(4>Mn,V)) = 0, (d in >)4.(u/)} = *r(" - u/), (12) 

where in the relevant limit T — > oo, St(0) = 1 and <5t(w) — > for u ^ 0. Phase noise modifies the corresponding 
relations for d; n i, 

(4l(^)^nl(^')) = l«inl| 2 , (JinlMJUw')) = - w') + (0(w)0(-o/)) |«i„l| 2 , (13) 

and adds correlations 

(di»i(w)dini (-"')) = - (£(w)£(-a/))a? n i- (14) 

We will parameterize the laser noise by an effective linewidth 7l(u>) given by the two-sided spectral density of frequency 
fluctuations, 7z,(w) = oj 2 ^/3(w)/3(— w)^. (For a laser with Lorcntzian lincshape, 7l is independent of frequency and 
represents the full width [7].) 



C. Transmission Spectrum 

(2) 

We now proceed to calculate the two-sided spectrum S R (u) of cavity transmission fluctuations and relate this to 
the spectrum of the mechanical oscillator's motion. The rate R at which photons are transmitted from the cavity is 
given in terms of the output field operator a out = aj n2 + - J~%a as R = a* ut a out . The fluctuations of this rate about its 
mean value R = a 2 ut are given by al ut a out - a„ ut rts (a* ut d ou t + aS out ), where d out = a out - a out and we are working to 
lowest order in d ou t/a ut- We assume, without loss of generality, that a is real. Defining e(u>) = d ou t(u>) +d out (—uj), we 
then have S R (uj)/R = (e(cj)e(- lu)). Using Eq. 7 to evaluate J ou t = d- m 2 + \J n/2d, we find e{uj) = e opt (uj) + e m och(^), 
where 

foptM = - (/t/2) x«(w)dini(w) + xl{-u)d\ nl (-u:) 

+ [1 - (k/2)x„(w)] <L2(w) + [1 - («/2) X ;(-w)] 4 2 (-^) (15) 

and 

emech(w) = ia-\/ K/2n(w) [c(oj) + c f (-w)] . (16) 

Here, e opt contains the intensity fluctuations due to photon shot noise or technical noise of the drive light, whereas 
£mech describes fluctuations in transmission due to atom-induced shifts of the cavity resonance. 

Using Eqs. 15 and 16, we can express the fractional fluctuations in transmitted intensity as the sum of three 
terms describing, respectively, the intrinsic optical fluctuations; the motion-induced fluctuations; and the correlations 
between the first two: 

S ( V(u)/Tf = sV(w) + + S^M- (17) 

Here, 

= (e opt He opt (-c))/i?= l/i? + 7L H|HH| 2 (18) 

represents the optical fluctuations that would be present even in the absence of optomechanical coupling, namely 
photon shot noise and laser phase noise (converted into intensity noise by the cavity). S^ ch (cj) is directly related to 
the spectrum S^(w) = (^X(lj)X(— lu)^ of the mechanical motion (evaluated below in Eq. 21) by 

SmechM = (emech(w)e mech (-w)) /R = Q 2 |IT(^)| 2 s£>(w). (19) 
Finally, the correlations between the mechanical motion and the optical noise are described by 

•Sfb M = (eoptMe mC ch(-w) + e me ch(w)e op t(-w)) /R 



( x ;H + 2 7L H|nH| 2 i?/«) 



(20) 



D. Displacement Spectrum 

We evaluate the two-sided spectrum of the mechanical oscillator's displacement (lu) = ( X(lu)X(—lu) ) using Eq. 



9. Applying the simplest possible model for the bath, namely quantum white noise [8] with ^c[ n (w)ci n (a;) ) = (ribath), 
we obtain 



7n 



(("-bath) + 1) IXmH-^)! + ("bath) I Xm* M I + 4 « ^tO-X^)^ + 8 [y L (w) R/ k] |w t all(w)| 



, (21) 



The first pair of terms (in square brackets) describes motion arising from the oscillator's coupling to the bath. The 
middle term describes motion induced by photon shot noise of the probe light, while the last term describes motion 
induced by laser frequency noise. 

Note that in the absence of optomechanical coupling (a = 0), the oscillator spectrum S x (uj) reduces to a pair of 
Lorentzians centered about ±cj t 



(("bath) + 1) |XmM| 2 + ("bath) |Xm(-^)| 2 



Q=0 



whose area is set by the bath temperature 



30 duj S { x\uj) 
co2tt XI 



2 (ribath) + 1. 



(22) 



(23) 



E. Comparison with Measured Spectra 



The spectra we measure are one-sided spectra, which we shall denote generically in terms of the two-sided spectra 

S^{co) by S(u) = S^(uj) + S^(-uj). The actual noise measured at the photodetector, normalized to the average 
2 2 

photocurrent I, is Si(w)/I = Sr(u)/R + Sdot- Here, Sdot accounts for imperfect quantum efficiency Q = 0.5(1), 
an excess noise factor F = 4.5(5) of the avalanche photodiode, and dark (primarily Johnson) noise, and is given by 



F-Q 2 Sj 
dct ~ Q R {QW 



(24) 



where Sj = 1.5(3) x 10 9 / (s 2 Hz) expresses the measured dark noise in units of equivalent photon rate. The spectra in 

Fig. 4 are taken with a photon rate R = 1.2 x 10 9 /s at the output of the cavity, which yields Sdot = 1-6(1) x 10~ 8 /Hz. 

2 

We can now write the measured spectrum as Si(u>)/I = S , mec h(w) + Sb g (w), where 



S hg (u) = S opt (uj) + Sfb(w) + Sdct- 



(25) 



From transmission noise S , op t(w) measured at large photon rate in the absence of atoms, we have determined the 
effective laser linewidth to be jL{^t) ~ 2n x 0.8(2) kHz at the trap frequency co t ~ 2tt x 500 kHz (well within our 3 
MHz lock bandwidth). At the photon rate R = 1.2 x 10 9 /s used in the spectra of Fig. 4, the associated phase- noise- 
induced intensity fluctuations are a factor of 1.7(5) below the photon shot noise level; correspondingly, they induce 
motion (included in our analysis) that is smaller than the zero-point fluctuations Xq. 

In fitting the measured spectra, we constrain F, Q, jl and R; we leave the dominant background noise contribution 
Sj free and obtain values consistent with the independently measured dark noise. 



F. Determination of Collective Temperature 

2 

Subtracting the background level 5b g (w) from the measured spectrum Si(u>)/I allows us to determine Sx((*>) from 
Eq. 21 and integrate it to find the occupation of the collective mode: 

(n)+l/2 = - ^ —S x (u)/X 2 . (26) 
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To determine (n) in Fig. 2, we use not the spectrum itself but the fractional variance o 2 = {I — I) 2 /I of the 

2 

measured photocurrcnt I oc R in a bandwidth B >• w t , which is related to the spectrum Si(lj)/I by 

, f B dui , , ,—2 f 00 do; _ . . f B dui „ , . ,„„. 

ff/ = y ~ y + y o ^^m. (27) 

We make in Eq. 19 for S , mcc h(w) the approximation H(uj) « II (wt) = (2/k)(>C+ — £-), which yields 

SmcchM « (2CJ/*;) 2 |£+ -£_| 2 5xH. (28) 

Integrating Eq. 28 allows us to obtain J" ^Sx(w) from erf using Eq. 27. The background noise term J Q S f^Sbg(w) 
on the right-hand side of Eq. 27 is independent of the occupation of the collective mode and is well approximated 
for the data in Fig. 2 by the variance erf oq measured in the long-time limit. In particular, we can find the change 
(n) — (n) in collective mode occupation between two different measurements erf, erf eq of the fractional transmission 
variance at fixed background noise by combining Eqs. 26-28: 

^ = 8(W«) 2 |£+ - £-| 2 (<n> - (n) cq ) . (29) 

The equilibrium occupation (n) < 3 of the collective mode C is small compared to the values (n) plotted in Fig. 2, 
where mode C is initially excited. Therefore neglecting (n) <C (n), and reexpressing Eq. 29 in terms of transmission 

2 

rate variances a 2 = (R — R) 2 /R = erf — -BSdct and cr^ g = erf cq — BXdct, we obtain 

^ - °l g = 8(GX / K ) 2 |£+ - £_ | 2 (n) . (30) 

III. ATOM NUMBER CALIBRATION 

We measure atom number [9] via the cavity shift Sun — NCQ, where C = N^ 1 J2iLi sm 2 (k^i). Allowing for the 
small but non-zero radial cloud size er r — 7(1) /an <C w, N represents an effective number of on-axis atoms. The 
cloud is long (w 1 mm) compared to the 5-/zm beat length between trap and probe, so that C = 1/2 in the absence 
of probe light. Displacement of the atoms by the probe light reduces C by at most 12% in our experiments, and we 
account for this effect. 



IV. THERMODYNAMIC TEMPERATURE 

The thermodynamic axial temperature, given by the mean single-atom vibrational occupation number (rii), is of 
interest for comparison with the bath temperature inferred from the fits in Fig. 4. We estimate (rii), in an ensemble 
of N = 1000(100) atoms, by ramping off the trap over 20 fis » 1/uJt while increasing the lattice depth of the probe 
to Up « 90hoj t . In the probe lattice, blue-detuned by A = +280 MHz from atomic resonance, axially cold atoms 
localize at positions x[ near the nodes. Before the cloud has time to expand radially, we determine the atom-probe 
coupling C = (sin 2 (fca^)) via the cavity shift, normalized by the shift measured beforehand in the 851-nm lattice at 
C = 1/2. From C w ((n^) + l/2)y/E T /U p we determine the mean vibrational level (n^(0)) = 6(2) in the final probe 
lattice. Note that this measurement provides only an upper bound on (rij) < (n'j) if the transfer into the deep probe 
lattice is not entirely adiabatic. 
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